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D^. Abstract 

bJQl We model the behaviour of a relativistic spherically symmetric shearing fluid undergo- 

ing gravitational collapse with heat flux. It is demonstrated that the governing equation 
for the gravitational behaviour is a Riccati equation. We show that the Riccati equa- 
tion admits two classes of new solutions in closed form. We regain particular models, 
obtained in previous investigations, as special cases. A significant feature of our solutions 
is the general spatial dependence in the metric functions which allows for a wider study 
of the physical features of the model, such as the behaviour of the causal temperature in 



^0 . inhomogeneous spacetimes. 

> 

X 1 Introduction 

An important application of Einstein's general relativity theory in relativistic astrophysics 
is the gravitational collapse of a radiating star. Since the first idealized model of a static 
spherical dust ball proposed by Oppenheimer and Snyder [Ij, many attempts have been made 
to describe more realistic situations. The derivation of the junction conditions for a radiating 
star with the exterior Vaidya metric was first obtained by Santos [2], and this formed the basis 
on which simple exact radiating models are constructed. The junction conditions were later 
generalised by Chan et al [5] to incorporate pressure anisotropy, and by Maharaj and Govender 
amongst others, to include the electromagnetic field. Particular solutions of the Einstein 
equations and boundary conditions have been found to describe the physical scenarios described 
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above. These solutions have been used to study the cosmic censorship hypothesis, and various 
physical features including the adiabatic stability of the model, surface luminosity, and the role 
of relaxational effects on the thermal evolution in the causal thermodynamic theory. 

In the absence of general techniques to analyse the nonlinear coupled partial differential 
equations in Einstein's theory, many formulations by researchers have sought to simplify the 
system of equations by considering the special case of shearfree collapse. A successful attempt 
at an exact model was made by Kolassis et a/ [S], where the assumption of geodesic motion of 
the fluid particles led to a considerable simplification of the field equations. Herrera et al [6j| 
solved Einstein's equations, and reduced the junction condition at the boundary, for shearfree 
collapse, to a nonlinear ordinary differential equation by requiing that all Weyl tensor com- 
ponents vanish. It was initially thought that this model could only be studied approximately. 
However this conformally fiat model was later solved exactly by Maharaj and Govender [7j by 
introducing a transformation that linearised the boundary condition. Herrera et al [8] subse- 
quently generalised the Maharaj and Govender [7] model and showed that other solutions to 
the linearised boundary condition existed. In a recent treatment, Misthry et al [9j extended the 
conformally fiat model to an exact nonlinear regime by transforming the governing equation to 
an Abel equation of the first kind. These solutions may be used to study realistic behaviour of 
the gravitating star with nonlinear boundary conditions. 

A natural extension of shearfree models is to include the effects of shear and pressure 
anisotropy. The origins and effects of anisotropy were first investigated by Herrera and Santos 
[To] and later by Chan et al [11] and Herrera et al [12]. The formulation of a shearing model 
with pressure anisotropy leads to a nonlinear partial differential equation at the boundary; 
this presents a formidable difficulty in obtaining solutions in closed form. Earlier treatments 
have been qualitative in nature without exact analytical solutions, and hence researchers have 
used numerical techniques to study the physical behaviour of the model. Noguiera and Chan 
|13j . assuming separable forms for the metric functions, reduced the boundary condition to 
a nonlinear ordinary differential equation. They then investigated the equation numerically 
and a detailed study of the physical features of the model was performed. Recently Naidu et 
al [Hj obtained the first exact analytical model with nonzero shear by considering geodesic 
motion of the fluid particles. Their particular solution, however, has singularities at the stellar 
core. Maharaj and Misthry |15] obtained two classes of exact solutions, nonsingular at the 
centre, in which the Naidu et al [H] model is contained as a special case. Note that the spatial 
components of the metric functions have been restricted to particular forms in the Maharaj 
and Misthry models. This restriction limits the investigation of the relativistic effects such 
as the relaxation time scales on the model, and consequently it is desirable to obtain a wider 
class of solutions. 

The main objective of this paper is to carry out a systematic study of the governing equation. 
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at the boundary, for the shearing collapse of a compact radiating stellar fluid model. We seek 
to obtain a general class of nonsingular exact solutions which allows for flexibility in the choice 
of physical parameters required to investigate the physical features of the model. It is not 
desirable to eliminate the inherent nonlinearity at the boundary; instead we seek to transform 
the governing equation to a familiar form, namely the Riccati equation. In Section 2, we 
formulate the model using the Einstein field equations together with the junction conditions. 
In Section 3, we first transform the governing equation into a Riccati equation, and then propose 
two transformations which lead to separable equations. Two new classes of exact solutions are 
found. It is shown that solutions found earlier with nonzero shear are contained in our models. 
In Section 4, we integrate the truncated form of the Maxwell-Cattaneo heat transport equation 
and obtain an explicit form for the causal temperature. We present profiles for the causal and 
acausal temperatures and briefly discuss the features of the graphs. 



2 Formulation of the model 

We seek to model a spherically symmetric star undergoing radiative gravitational collapse with 
nonzero shear in the context of general relativity. The line element describing the gravitational 
field for the interior spacetime is taken to be 

ds"^ = -dt^ + B'^dr'^ + Y^de"^ + sin^ ed(f)'^) (1) 

where B and Y are functions of both the temporal coordinate t and radial coordinate r. For 
this model, the fluid four- velocity = 6q is comoving. The fluid four- acceleration vector ii"', 
the expansion scalar 0, and the magnitude of the shear scalar a, respectively, are given by 

ii" = (2a) 

e = f I + 2^^ (2b) 

(2c) 

for the line element ([T]). We observe that the particle trajectories of the collapsing fluid are 
geodesies because ii"- = 0. However note from ([21) that the fluid expansion and the shear may 
be nonzero in general. A similar analysis was performed by Kolassis et al [5] when the shear is 
vanishing; in that case it was possible to solve the boundary condition and the field equations 
and obtain an exact solution. The Kolassis et al |5j model has a Friedmann limit and the dust 
cosmological model is regained. It would be interesting to compare the temporal evolution of 
the model when shear is present. Investigation of the behaviour of the temperature in causal 
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thermodynamics for geodesic motion has revealed higher central temperatures than the Eckart 
theory as established by Govender et al [H]. In this study we seek to incorporate the effects of 
shear in the model. 

The energy momentum tensor for the interior spacetime has the form 

Tab = {p+ P)UaUb + PQab + QaUb + QbUa + T^ab (3) 

where p is the density of the gravitating fluid, p is the isotropic pressure, qa is the heat flux 
vector and i^ab is the stress tensor. These quantities are measured relative to the four-velocity 
u. The stress tensor can be written explicitly as 

VTab = {Pr - Pt) (naUb - ^^afe j (4) 



where Pr is the radial pressure, pt is the tangential pressure, and n is the unit radial vector 



orthogonal to u. Hence we have = 45". The isotropic pressure 



P=l{Pr + 2pt) (5) 

relates the radial pressure and the tangential pressure. 

It is possible to write the Einstein field equations as the set 

BY 1 Y^ I f Y" Y'^ B'Y'\ , , 

P = ^BY^Y-^^Y-^-T^i'Y+Y^-'BY) ^''^ 

Y y\ 1 1 

Pr = |-V"l^ + Vj+5^(,l^J"l^ ^^^^ 

B'Y' Y"\ _ _ 




B^ \ BY Y 
2 f Y' BY , 

for the spacetime ([I]) and the matter distribution ([3]). The fluid pressure is anisotropic and 
the heat flux q°- = (0, q, 0, 0) has only a nonvanishing radial component. We observe that if 
functional forms for the gravitational potentials B and Y are given then expressions for the 
matter variables p,pr,pt and q immediately follow from ([6]). 

The exterior spacetime, describing the region outside the stellar boundary, is described by 
the Vaidya metric 

ds^ = -(l- ^^"j dv^ - 2dvdR + R\de^ + sin^ ed(j)^) (7) 
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where m{v) denotes the mass of the fluid as measured by an observer at infinity. The metric ([7]) 
describes coherent null radiation which is fiowing in the radial direction relative to the hyper- 
surface E which represents the boundary of the star. The matching of the exterior spacetime 
with the interior spacetime leads to the following set of junction conditions on the hypersurface 
E: 

1/2 



dt 

{Prh 



1 - 

!(■ 

{qBh 



2m ^ JR^ 
Rs dv 



dv 



J E 



(8a) 
(8b) 
(8c) 
(8d) 



The nonvanishing of the radial pressure pr at E leads to an additional differential equation, 
namely the boundary condition (l8dl) . which has to be satisfied together with the field equations 
This condition was first established by Santos [2] for shearfree spacetimes, and extended to 
spacetimes with nonzero shear by Glass [T7] and Maharaj and Govender [IB], amongst others. In 
a recent investigation by Di Prisco et al [T9] , the matching conditions applicable to spherically 
symmetric gravitational collapse with dissipation and nonzero shear have been generalised to 
include nonadiabatic charged fiuids. 



3 Solution of the governing equation 

The junction condition (prOs = {.(iB)t. becomes 

2YY + Y^-Xl + lyf - 2^YY' + 1 = (9) 

which follows from ([6]). Equation governs the gravitational behaviour of a radiating star with 
anisotropic pressure and nonzero shear. To complete the description of the gravitational be- 
haviour of the model we need to integrate the junction condition (Q; this will lead to functional 
forms for the metric functions B{r,t) and Y{r,t). Exact solutions for the junction condition 
([9]) have been extremely difficult to obtain due to the nonlinear nature of the equation. In an 
earlier study of a shearing radiating model, Noguiera and Chan [13] used numerical techniques 
to obtain approximate solutions. Ideally an exact solution is desirable in terms of elementary 
or special functions. An exact solution for this physical model was obtained by Naidu et al 
[14] in terms of the elementary functions. This class of solution is singular at the stellar cen- 
tre. Maharaj and Misthry [15] extended the Naidu et al [H] model and showed that a wider 
category of solutions are possible; the singularities at the centre were shown to be avoidable. 
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We seek, in this treatment, to obtain a general class of nonsingular solutions which will 
allow for an investigation of the physical features of the model. Previous treatments were ad 
hoc. Our objective is to write ([9]) in a generic form, and then obtain solutions systematically. 
It is shown that the model leads to the formation of a Riccati equation where the potential B 
is the dependent variable. We present, in the following, a method of solving ([9]) exactly, and 
find several classes of solutions depending on the form of Y used. We rewrite iQ in the form 



B=\t7: + 7^ + 7^]B' + tt:B-— (10) 




which is a Riccati equation in the potential B. We demonstrate that two classes of solutions 
can be found for this Riccati equation. 

3.1 The first solution 

We seek solutions where F is a separable function of the form 

Y = R{r){t + af/'' (11) 

so that the temporal evolution of the model is specified. It is convenient at this point to 
introduce the transformation 

B = Z{t + a)2/3 (12) 

Then ( fTOj) can be written in the form 

(t + afl^Z = -^iZ^ - R'^) (13) 

Equation ( |T3i) is simple and separable with solution 

y^j., / l + ^(r)exp[3(t + a)VVi^] \ 
\l-g{r)e^Y>W + <^Y'VR]) 

where g{r) is related to an arbitrary constant of integration. Consequently the potential B can 
be obtained explicitly in the form 

ll-,(r)exp[3(t + a)V3/i?]J(^ + «) (^5) 



From f[TT|) and flTSl) we may write the interior metric ([T]) as 

g{r) exp[3(t + a)i/3/i?]y 



ds" = -de + (t + a)^/3 



R ^-77777 -I/// J Ctr^ + ^ + sm^ 

V 1 - ofr)exp 3ft + a)i/3/i? y ^ 



(16) 



which describes the interior spacetime of the radiating star. 
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The matter variables for the model are given by 



4 

3^2 



1 + 



3t^ 



gexp[3p/^/R] 



2/3(1 + ^exp[3tV3/i?]) 



8gexp[3t^^^/R] 



Pr 



Pt 



3t5/3i?(l _^2gxp2[3tl/3/i?]) 

4exp[3p/Vi?](l - gexp[3i^/^/R]) 
RR'i{l + g[exp 3tV3/i?])3 

4^exp[3p/Vi?] 

'i^/3R^{l + gexp[3i^/yR]^ 

2gexp[3i'/yR] T 
"i?2t4/3(i_^2exp2[3ti/3//?]) [ I 

2exp[3p/Vi?](l - ^exp [3i^/^/R]) 
RR'i{l + gexp [3t i/3//?])3 

4^exp[3p/Vi?](l - ^exp [3p/Vi?]) 
R'^RPil + ^exp [3t i/3/i?])3 



tV3 



3^ 

i?2 



(17a) 
(17b) 



2c/exp[3p/Vi?] m 
l-gexp[3W^/R]) ^ sTV^ 

p/3 i?2 



(17c) 
(17d) 



which satisfies the Einstein equations ([6]). For simplicity we have set t = t + a. 

We have obtained an exact solution to the Einstein field equations ([6]) that satisfies the 
boundary condition ([9]) for a radiating relativistic star. This is a new class of exact solutions 
where the spatial dependence in the function R is arbitrary. Consequently particular solutions 
found in the past can be shown to be contained in this class. We observe that when R = r + b 
the metric f|T6l) becomes 



-dt^ + {t + 



N4/3 



+ {r + b)\de^ + sin'^ 



' ^(r) exp[3(t + a)^/V(^ + ^)] 
^(r)exp[3(t + a)V3/(r + 6)] 



dr' 



(18) 



Thus we have regained the Maharaj and Misthry [T5] model which is regular at the stellar 
origin. The Naidu et al [Hj model is regained from (fT6|) when a = and 6 = 0. Other choices 
for the function R are clearly possible: the choice should be such that the model remains regular 
at the centre and the model is well behaved in the interior. We further observe that the model 
yields the Friedmann dust model when g = 0. In this case we can find coordinates such that 



ds^ = -de + ^"'[dr^ + r\de^ + sin^ Od^p")] 



(19) 



for which the heat fiux vector vanishes and Pr = Pt = ^ with p = ^ in geodesic motion. 
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3.2 The second solution 



Other solutions to the Riccati equation ( JTOl) exist in closed form but these are difficult to find 
in practice. It is possible to find a second class of solutions to ffTOl) by assuming that 



Y = R{r)(t + a) 



In this case we introduce the transformation 



B = (t + a)Z 



Then f l20|) can be written in the form 



{t + a)Z 



R'' 



2RR' V + i 

Equation (l22l) is a simple equation, separable in the variables Z and t, with solution 



(20) 



(21) 



(22) 



R' 



VW+i 



l + hir)it + ay 



\ 



(23) 



\l-h{r){t + a)^ « / 



where h{r) is an arbitrary function of integration. Therefore we can express the metric potential 
in the form 

R' l + h{r){t + ay^ ^ 

Jd — 



Then from (l20i) and we obtain the interior metric 



(t + a) 



(24) 



-dt^ + {t + af 



i?2 + 1 

which describes the stellar interior 



R 



12 



( 



1 + /i(r)(t + a 



2^ 



rfr' + i?"(c/r + sin 



(25) 
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The Einstein field equations ([6]) then imply the matter variables which are given by 

-\ 2 



P = 7^ 



2 



1 + 



- (3i?2 + 1) 



1 - hi 



Rt\ 1 + ht 



2h^WTl t 



RR' 



t2 \l-h^t 



,/ -R^ + l 

?2 , in2 iV~r2- 



+ 



4(/?2 + i)^ t 



h' hint 



i?2 + 1 i?9 



(26a) 



Pr 



Pt 



1 R^ + 1 
P R^ 



\-hV ^ 



- 1 



\\\hV « / 
2h^/WTTi 2\/^ 



i?t2 1 _ /^2^ 2V^ 



1 + 



+ 1 _ 1 

RR' R 



1 



1 

12 



1 -ht 



\l + hV^ J 



Ah^{R^ + l)t 



I R'^ + 1 



'fl2 + l 



i?2i?'t2 1 _ /i2^ 2V^ 1 - /it 



2(i?2 + i)t ■ 



RR'P [l + hi 



h' 



h{R^ + l)\ni 



R9 I l^hp^ 



I R^ + 1 



4/l(i?2 + l)3/2^~V {l-hV ^ 



/ r'^+1 



R^R'"^ \ l + hV^ 



(26b) 



(26c) 



(26d) 



where we have set t = t + a for convenience. 

We have generated another exact solution to the field equations ([6]) that satisfies the bound- 
ary condition ([9]). Again the spatial dependence in the function is arbitrary and the model is 
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regular at the centre. We observe that the solution contains an exponential temporal depen- 
dence on the spatial function R{r) of the form tV^R^^^/^. Such solutions are difficult to 
interpret but may be relevant in describing new physical models in the strong gravity regime 
for gravitational collapse. 



4 Causal temperature 

The simple forms of the solutions found in this paper, in particular the first solution, permit a 
study of the physical features. We consider briefiy the relativistic effects on the temperature. 
For a shearing superdense matter distribution, we employ the Maxwell-Cattaneo heat transport 
equation to investigate the causal thermodynamical behaviour of the model. In the absence of 
rotation and viscous stress this is given by the truncated version 

rh% + qa = -K{hybT + Tua) (27) 

where r is the relaxation time, n is the thermal conductivity, and hab = Qab + UaUb projects into 
the comoving rest space. When r = we regain the acausal Fourier heat transport equation. 
For our model equation (|27|) may be written as 

T = -- I {qB)'Bdr - - [ qB^dr (28) 



K J K ^ 

describing the evolution of the causal temperature. It has been shown in previous investi- 
gations of relativistic stellar models that the effect of the relaxation time r, on the thermal 
evolution, plays a significant role in the latter stages of collapse [201 ED [221 [23|- Naidu et al 
[H] showed that in the presence of shear stress, the relaxation time decreases as the collapse 
proceeds and the central temperature increases. The particular form of the relaxational time 
r is dependent on the physical constraints of the model during the latter phases of collapse. 
We observe that since our solutions have elementary functions with an arbitrary form for the 
spatial component, it is possible to integrate ( 128|) for different choices of r. In particular, the 
effect of decreasing relaxation time with decreasing radius and higher central temperature is 
possible by incorporating a varying function r. 

In this study we set r and k to be constant in order to examine the causal and acausal 
behaviour in the first solution (ITS]) . We need to choose particular forms for the arbitrary 
function R{r) to complete the integration. As a first example we set i? = r + 6 to obtain for 
the temperature 



At 



9Kt 



2 



Log 



[l+gexp[3t^^^/r + b]y 
-l + gexp[3i^/^/r + b]y_ 



3Kt5/3 \{r + b){l + gexp[3ty^/r + b])J 3Kt 



=Tanh-^(^exp[3t '/7r + b]) + f{§9) 
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where we have kept g constant, and f{t) is a constant of integration related to luminosity as 
observed by a distant observer. As a second example we set R = to obtain the temperature 



4r 

T = -^(Log 



;i + ^exp[3t^/Ve''])= 



4r ( ^exp[3F/3/e'^] ^ 4 _ , -1/3 



.Tanh-^(^exp[3t ^/7e"]) + f{t) (30) 



3Kt 5/3 \^ e^' ( 1 + ^ exp [3t ^3 /g^] ) y 3^^ 

Consequently it is possible to find analytic forms for the causal temperature in terms of elemen- 
tary functions as shown in fl29|) and fl30|) . We regain the noncausal (Eckart) temperature profiles 
when r = 0. Our simple forms for T assist in studying the evolution of a radiating star in dif- 
ferent time intervals. These models provide examples of temperatures where inhomogeneity is 
directly related to dissipation. 

It is possible to qualitatively distinguish the causal and acausal temperatures for the region 
between the centre and the surface of the star. In Figure [H we provide plots of the causal 
(solid line) and Eckart (broken line) temperatures against the radial coordinate on the interval 
< r < 1, where we have selected r = 1 for simplicity. This figure has been generated with the 
help of Mathematica. We observe that the temperature is a monotonically decreasing function 
as we approach the boundary from the stellar centre. Also, it is immediately clear that the 
causal temperature is everywhere greater than the acausal temperature in the interior of the 
star. At the boundary E, however 

T{t,rY)causal = T{t 

) ^T,) acausal 

(31) 

This simple figure has been generated by assuming a particular constant value for the relaxation 
time T and the thermal conductivity k. Changing the magnitude of these values would produce 
a change in the separation of the curves but the results do not change qualitatively. For example, 
in Figure El we provide a plot of the causal (solid line) and Eckart (broken line) temperatures 
for r > 0. We note that in this case both temperatures decrease more rapidly as we approach 
the boundary; the value of r affects the gradient of the temperature. As indicated previously 
it is possible for the relaxation time r to vary. The choice for r should be dictated on physical 
grounds, e.g. rate of particle production at the stellar surface. 
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